Violation of Ioffe-Regel condition but saturation of resistivity of the high T c cuprates 
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We demonstrate that the resistivity data of a number of high T c cuprates, in particular 
La2- a: Sr a; Cu04, are consistent with resistivity saturation, although the Ioffe-Regel condition is 
strongly violated. By using the f-sum rule together with calculations of the kinetic energy in the 
t — J model, we show that the saturation resistivity is unusually large. This is related to the strong 
reduction of the kinetic energy due to strong correlation effects. The fulfilment of the Ioffe-Regel 
condition for conventional transition metal compounds is found to be somewhat accidental. 



Some strongly correlated metals [1-4], in particular 
many high T c cuprates [5-14], show an exceptionally 
large resistivity, p, at high temperatures T. Thus p can 
reach values of several mi7cm. This is in strong con- 
trast to almost all other metals. Typically, a metal with 
a very large resistivity shows resistivity saturation [15]. 
Thus when p reaches values of the order 0.1 mficm, the 
slope of p(T) is typically reduced substantially, although 
not necessarily to zero. This happens when the apparent 
mean free path I becomes comparable to the separation, 
d, of two atoms, the Ioffe-Regel condition [16]. This kind 
of behavior has been observed for many metals and it 
used to be considered a universal behavior [17]. 

The metals mentioned above and the alkali-doped 
fullerenes [18] are apparent exceptions to this saturation 
behavior, with the cuprates forming a particular chal- 
lenge [19,20]. The deviation from the "universal" behav- 
ior can be illustrated by considering Ioffe-Regel condition 
for La2_ a; Sr a ;Cu04. We assume a "large" Fermi surface 
of cylindrical shape, containing 1 — x electrons (or 1 + x 
holes) [21]. Assuming that I = a, where a is the lattice 
parameter in the CuC>2 plane, we obtain the saturation 
resistivity 

0.7 

p sat = mOcm. (1) 

VI ± x 

If we instead assume a "small" Fermi surface, y/l ± x in 
the denominator is replaced by y/x. In cither case, the 
experimental resistivity for small x is much larger than 
the saturation resistivity above, supporting the conclu- 
sion that there is no saturation in these systems. 

Here we demonstrate that the experimental data for 
the strongly correlated high T c cuprates are, neverthe- 
less, consistent with saturation. The saturation resis- 
tivity can, however, be much larger than the Ioffe-Regel 
value. This is not entirely surprising, since the Ioffe-Regel 
condition is based on a semiclassical picture, which is not 
valid when I ~ d. The good agreement between the Ioffe- 
Regel condition and the saturation for weakly correlated 
transition metal compounds is therefore somewhat acci- 
dental. The large saturation resistivity is shown to be due 
to a large reduction of the kinetic energy, due to strong 
correlation effects. Since the alkali-doped fullerenes ap- 
pear to lack saturation [22], this puts these compounds 
in a special class, different from the cuprates. 



To discuss resistivity saturation, we use the f-sum rule 
in a form appropriate for the models discussed here (only 
nearest neighbor hopping and no on-site matrix elements 
of the current operator) [23,24] 



2 r < w 1 

- / a{uj)duj = -- 

7T Jo 2 



(Tk), 



(2) 



where <j(ui) is the optical conductivity, d is the separa- 
tion of the sites, f2 is the volume of a unit cell, N is the 
number of unit cells and (Tk) is the expectation value of 
the kinetic energy operator. The prefactor 1/2 refers to 
the two-dimensional case and is replaced by 1/3 for the 
three-dimensional case. To obtain an approximate up- 
per limit to the resistivity, we assume that the (Drudc) 
peak at uj = has been smeared out and that a(uj) is a 
smooth function. The removal of the ui = peak may be 
due to any scattering mechanism, electron-phonon [20], 
electron-electron or impurity scattering and, depending 
on the system, it may or may not happen at large val- 
ues of T. We emphasize that our theory below does not 
depend on any particular scattering mechanism. We fur- 
thermore assume that <r(u) = for \ui\ > W, where W 
is the one-particle band width. This should be a good 
approximation, since \u>\ > W would involve multiple 
electron- hole pair excitations and have a small weight. 
Explicit calculations for a model of transition metal com- 
pounds and for the t — J model confirm this. If a(u) had 
a box shape, cr{uj) = u(0) for < \u>\ < W, the integral 
over cr(u>) would be a(0)W and we could write 



(t(oj = 0) = — — f°° u 
K ' Wj 



(u))du>, 



(3) 



with 7=1. For a more realistic shape of (t(u>), where 
u{ui) has a maximum for u = 0, it follows that 7 > 1. 
Following the considerations for a model of transition 
metal compounds with electron-phonon scattering [25], 
we may set 7 = 1.9. In the t — J model, considering 
electron-electron scattering, we find that <j(u>) is more 
concentrated to small u>, giving a larger value of 7. If 
the electrons would tend to localize, 7 could be much 
smaller than unity. This should, however, not happen 
at the large values of T considered here. Measurements 
[26] of cr(u>) for La2- x Sr x Cu04 at large T show a similar 
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behavior as assumed in Eq. (3), apart from some sharp 
phonon structures at finite u. 

We first consider noninteracting electrons in a band 
with the orbital degeneracy n. The scattering is assumed 
to be due to the electron-phonon interaction. For T <C W 
and neglecting the influence of the electron-phonon cou- 
pling on the kinetic energy, 

(T K ) = 2n eN(e)de = -2naWN, (4) 

J-W/2 

where N(e) is the density of states per orbital and spin, 
\x is the chemical potential and a depends on the shape 
N(e) and the filling. For semi-elliptical N(e) and filling 
0.4, a = 0.10. This result is relatively independent of the 
specific shape of N(e), and using, e.g., a constant N(e) 
leads to a = 0.12. For a three-dimensional system we 
then obtain 

irn d 3 e 2 . . 

^ = i^y' (5) 

which provides an approximate upper limit to the resis- 
tivity for T <C W. Saturation happens if p(T) grows so 
rapidly that this approximate upper limit is approached 
for experimentally accessible values of T. This happens 
for some transition metal compounds, with the A15 com- 
pounds, such as NbsSb, being particularly striking ex- 
amples. The important orbitals in NbsSb are the Nb d- 
orbitals, leading to n = 5. Considering the A15 lattice of 
Nb 3 Sb, we obtain the saturation resistivity 0.14 mf2cm, 
which compares well with the experimental saturation 
resistivity of about 0.15 mffcm [15]. This illustrates the 
usefulness of this approach for obtaining a saturation re- 
sistivity. 

We next consider the strongly correlated cuprates. To 
calculate the kinetic energy (Tk), we need an appropriate 
model for these systems. The transport properties should 
mainly be determined by the Cu-0 antibonding band of 
Cu x 2 — y 2 and O 2p character. Thus we consider one 
orbital per site and the orbital degeneracy n — 1. These 
orbitals are put on a two-dimensional square lattice, and 
we consider a a layered system with the lattice param- 
eters (a,a,2c). We further assume that the Coulomb in- 
teraction U between two electrons on the same site is so 
large that double occupancy of a site can be neglected. 
This can be described by the t — J model [27]. The 
properties of this model have been studied extensively 
by Jaklic and Prelovsek [24]. The system has the hole 
doping x (> 0). We emphasize that the t — J model 
is only used to estimate (Tk)- We do not exclude other 
scattering mechanisms; we just neglect their contribution 
to the kinetic energy. 

For simplicity, we first consider J = 0. By expanding 
in 1/T, we obtain the large T limit of the kinetic energy 
as 

Tr T K e~ H / T 1 Tr T 2 
(Tk) = - 

Tr e~ H / T ~ T ^ Tr 1 doubWo occ ' 

(6) 



where the right hand site is evaluated for states without 
double occupancy. Since the expression contains T K , this 
involves a hole hopping to a neighboring site and back. 
In the t — J model this is only possible if the neighbor- 
ing site has no holes. The probability for having a hole 
or no hole on a given site is x and (1 — x), respectively. 
The probability for hopping between two nearest neigh- 
bor sites is then approximately x(l — x). Summing over 
the four nearest neighbors, this gives the large T limit 

(T K ) = -ANt 2 ^^!. (7) 




x 

FIG. 1. The kinetic energy Tk(x) as a function of doping 
x in the t — J model for J = 0. The crosses show the result of 
exact diagonalization for 4 x 4 two-dimensional cluster [30], 
the full curve Eq. (8) and the broken curve the result for spin- 
less fermions [28,29]. The figure illustrates the approximate 
x(l — x) behavior of the kinetic energy. 

Using similar arguments, we can obtain a rough esti- 
mate for the T = kinetic energy 

[(T k )]t=o = -4te(l - x)N. (8) 

This result is shown by the full curve in Fig. 1. A 
similar result, but with a smaller coefficient, can be ob- 
tained variationally, by assuming that the holes behave as 
spinlcss fermions in a ferromagnetic background (broken 
curve in Fig. 1) [28,29]. These results can be compared 
with results from an exact diagonalization [30] of a 4 x 4 
t — J model with J = (crosses in Fig. 1). The figure 
illustrates the approximate x(l — x) dependence of the ki- 
netic energy. In the following we use the t — J model with 
J/t = 0.3 and t = 0.4 cV [24]. We solve the model, using 
the finite T exact diagonalization technique of Jaklic and 
Prelovsek [31]. We consider a finite cluster of 4 x 4 atoms 
and use the grand canonical ensemble. The resulting nu- 
merical results for the kinetic energy for small values of 
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T and x can be approximately expressed in the form of 
Eq. (8), but with the prefactor 4 reduced to about 3.4. 
For the values of T of interest here, the T dependence 
of (Tk) is weak, and it is neglected in the following. If 
the system has a strong clcctron-phonon interaction, the 
kinetic energy should be somewhat further reduced. 

Combining Eqs. (2, 8), we predict that the right hand- 
side of Eq. (2) should be proportional to x, for small x. 
The effective number of carriers per Cu atom N e f / (x) can 
be estimated integrating the optical conductivity up to 
cj c = 1.2 eV, as suggested by Yamada et al. [32]. This also 
includes some contribution from interband transitions, 
and Yamada et al. [32] therefore subtracted N e ff(0), 
which should contain only interband transitions. They 
found that N e ff(x) — N e fj(0) is indeed approximately 
proportional to x, up to x = 0.12. Based on the absolute 
numbers by Uchida et al. [33], we find that our factor 
of proportionality in Eqs. (2, 8) agrees with experiment 
to within about 20-30 %. Due to the problems of re- 
moving the interband transitions, this should be a good 
agreement. 

We consider conduction in the ab-plane of a cuprate 
and assume that the u> = peak is gone. Insertion of Eq. 
(8), but with the prefactor 3.4, in Eqs. (2, 3) gives an 
approximate upper limit to the resistivity for small and 
intermediate values of T 



0.4 



x(l — x) 



mficm, 



(9) 



where we have used distance c = 6.4 A between the Cu02 
planes. Inserting the calculated kinetic energy for the 
t — J model would lead to a somewhat larger resistivity 
for small x and it would introduce a weak T dependence. 
The saturation resistivity in Eq. (9) is much larger than 
Eq. (1) for small x. It illustrates that one should ex- 
pect saturation at a substantially larger resistivity than 
predicted by the Ioffe-Regel condition. The saturation 
resistivity in Eq. (9) is also much larger than the result 
derived [25] for a model of weakly correlated transition 
metal compounds 



0.2d 0.5 



(10) 



where d is expressed in A. This is partly due to the de- 
generacy being just n = 1 for the t — J model but n = 5 
for our model of the weakly correlated transition metal 
compounds. Furthermore, the strong correlation drasti- 
cally reduces the kinetic energy in the cuprates, which 
shows up as the factor x(l — x) in the denominator of 
Eq. (9). Finally, the large value of the lattice parameter 
c also increases p for the cuprates. 

We now consider the T dependent resistivity for dif- 
ferent doping. The experimental results of Takagi et al. 
[6] are shown in Fig. 2. For T < 300 K, more recent 
results have been obtained by Ando et al. [7]. Their re- 
sistivities are qualitatively similar but generally smaller 



than the results in Fig. 2. The results by Takagi et al. [6] 
are compared with the saturation resistivity in Eq. (9). 
For all values of x, the experimental resistivity is smaller 
than the saturation resistivity predicted here. The ex- 
perimental results do therefore not demonstrate absence 
of saturation. Actually, for x = 0.04 and x = 0.07 the 
data show signs of saturation for values not much smaller 
than the expected saturation resistivities, but at much 
larger values than the saturation resistivity predicted by 
the Ioffe-Regel criterion (Eq. (1)) (indicated by an ar- 
row). This shows that the Ioffe-Regel criterion is in- 
valid for these systems. It is interesting to plot xp(T) 
[34]. In such a plot the saturation resistivity becomes 
almost independent of x. Indeed the curves of xp(T) for 
x = 0.04 and x = 0.07 fall almost on top of each other, 
suggesting that close to saturation there is a scaling of 
p(T) by 1/x, as implied by the arguments above. For 
x = 0.15 and x = 0.34 one would have to study much 
higher (and experimentally unaccessible) values of T to 
determine whether or not there is saturation. 
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FIG. 2. The resistivity as a function of T for 
La2-a;Sra;Cu04 [6] and the saturation resistivity according 
Eq. (9) for x = 0.04 (full curve) x = 0.07 (broken curve) 
x — 0.15 (dotted curve) and x — 0.34 (chain curve). The hor- 
izontal arrow shows the saturation resistivity expected from 
the Ioffe-Regel condition. The figure illustrates that there are 
signs of saturation for small x at roughly those resistivities 
where saturation is expected. For larger x, much larger T 
would have to be considered to test whether there is satura- 
tion. 

I^-aSrajCuO^ has a structural transition from an or- 
thorhombic to a tetragonal phase. This transition hap- 
pens, however, at a somewhat lower temperature [32] 
than the break in the slopes of the curves for x = 0.04 
and x = 0.07 in Fig. 2. 

We have also considered the resistivity 
of YBa 2 Cu 3 6+5 [8], Nd 2 _ x Ce x Cu0 4 -j, [9], 
Bi 2 Sr 2 Ca 1 _ x Y x Cu 2 8+!/ [10,11], Tl 2 Ba 2 Cu0 6+ ^ [12] 
and HgBa 2 Ca„_iCu„0 2 „ +2 +a; [13]. In all these cases 
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we find that the experimental resistivity does not ex- 
ceed the expected saturation resistivity. In the case of 
Bi2Sr2Cai_ x Y x Cu208+j/ saturation has been reported 
[11]. This provides further support for our theory. 

Above we have provided an upper limit to the resis- 
tivity based on the f-sum rule and the calculated kinetic 
energy and without specifying a scattering mechanism. 
We have also calculated the resistivity of the t — J model 
directly, although the results are not very accurate for 
the small clusters considered here. For x = 0.15, the 
calculated resistivity is of the right order of magnitude, 
while for x — 0.07 and x = 0.04 it is substantially smaller 
than the experimental results in Fig. 2. This suggests 
that there is an additional scattering mechanism, beyond 
the electron-electron scattering included explicitly in the 
t — J model, which becomes important for x = 0.04 and 
x = 0.07. This additional mechanism apparently makes 
the resistivity increase so rapidly for small T that the 
approximate upper limit (9) is approached and the re- 
sistivity shows sign of saturation. The estimate of the 
upper limit (9) is nevertheless correct, unless the addi- 
tional scattering mechanism influences the kinetic energy 
in a substantial way. 

Above we have discussed the situation when T is much 
smaller than the band width, which applies for values of 
T which can be reached experimentally. The kinetic en- 
ergy Tk then has a weak T dependence, leading to an 
upper limit for the resistivity which is only weakly T de- 
pendent. As can be seen from Eq. (7), however, the T 
dependence becomes strong for very large T (> W/A). 
Inserting Eq. (7) into Eqs. (2,3) then leads to an upper 
limit for p(T) which grows linearly with T . A similar 
result would be found even for noninteracting electrons 
scattered by impurities, i.e., for a T- independent mecha- 
nism, where the T dependence can then be traced back 
to the T dependence of the Fermi functions. This illus- 
trates that "saturation" does not necessarily imply that 
the resistivity becomes T independent, just that its slope 
is reduced substantially at large T [17]. The effects dis- 
cussed in this paragraph, however, refer to very large 
values of T, of the order of 10000 K. They are therefore 
only of academic interest, unless the band width is very 
small. 

To summarize, we have shown that the experimental 
resistivity data for a number of cuprates are consistent 
with resistivity saturation. The saturation resistivity is, 
however, much larger than for, e.g., A15 compounds, due 
to the strong suppression of the kinetic energy in the 
strongly correlated cuprates. The Ioffe-Regel condition 
is therefore violated for a many of these systems. 
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